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We consider a general relativistic zero-pressure irrotational cosmological medium perturbed to the 
third order. We assume a flat Friedmann background but include the cosmological constant. We 
ignore the rotational perturbation which decays in expanding phase. In our previous studies we 
discovered that, to the second-order perturbation, except for the gravitational wave contributions, 
the relativistic equations coincide exactly with the previously known Newtonian ones. Since the 
Newtonian second-order equations are fully nonlinear, any nonvanishing third and higher order 
terms in the relativistic analyses are supposed to be pure relativistic corrections. In this work we 
derive such correction terms appearing in the third order. Continuing our success in the second-order 
perturbations we take the comoving gauge. We discover that the third-order correction terms are 
of (p„-order higher than the second-order terms where ^„ is a gauge-invariant combination related 
to the three-space curvature perturbation in the comoving gauge; compared with the Newtonian 
potential we have (5<I> ~ to the linear order. Therefore, the pure general relativistic effects are 
of ip^-order higher than the Newtonian ones. The corrections terms are independent of the horizon 
scale and depend only on the linear order gravitational potential (curvature) perturbation strength. 
From the temperature anisotropy of cosmic microwave background we have 4^ ~ i5$ ~ ~ 
10"^. Therefore, our present result reinforces our previous important practical implication that near 
current era one can use the large-scale Newtonian numerical simulation more reliably even as the 
simulation scale approaches near (and goes beyond) the horizon. 



I. INTRODUCTION 

In our previous works [1,2] we have proved that in 
the zero-pressure irrotational cosmological medium the 
relativistic second-order scalar-type perturbation equa- 
tions coincide exactly with the known ones in Newto- 
nian theory. This result shows a continuation of the 
relativistic-Newtonian correspondences of the cosmolog- 
ical zero-pressure medium previously shown in the back- 
ground world model by Friedmann in 1922 [3] and by 
Milne in 1934 [4], and in the linear perturbation by Lif- 
shitz in 1946 [5] and by Bonnor in 1957 [6]. History shows 
that both for the background and for the linear perturba- 
tion the equations were first derived in Einstein's gravity 
[3,5] and later followed by studies in Newton's gravity 
[4,6]. In the case of second-order perturbations the New- 
tonian result was known first [7]. Since the Newtonian 
results are not supposed to be reliable as the scale ap- 
proaches the horizon, our result has a practical impor- 
tance by showing that even to the second-order pertur- 
bations we can use Newtonian equations in all scales. In 
this work we will extend the situation to the next order in 
relativistic perturbation which, as we will explain shortly, 
can be regarded as pure relativistic corrections. Even in 
the linear perturbations the presence of the tensor-type 
perturbation (gravitational waves) can be regarded as a 
pure relativistic effect. To the second order, additionally, 
the gravitational waves and the scalar-type perturbation 
are coupled, thus one can source the other. This can be 
regarded as another pure relativistic correction. 



In this work we will derive pure relativistic corrections 
of the scalar-type perturbation which appear in the third- 
order perturbations. The situation we have is that to the 
second order Einstein's theory gives exactly the same re- 
sult as the Newton's theory whereas in the Newtonian 
case such second-order equations are in fact exactly valid 
to fully nonlinear order. Therefore, any nonvanishing 
third and higher order perturbations in Einstein's grav- 
ity are supposed to be the pure relativistic corrections. 
We will present the scalar-type perturbation equations 
to the third order and the gravitational wave equation 
to the second order. This is enough because the previ- 
ously known gravitational wave equation will already get 
the correction terms even in the second order, whereas for 
the scalar-type perturbation new correction terms appear 
only in the third order compared with the previous stud- 
ies [2]. If we include the gravitational waves, the gravi- 
tational waves only to the second order will be needed to 
make the third-order scalar type perturbation equations 
complete. 

We will ignore the vector-type perturbation in this 
work, thus consider only irrotational perturbation, be- 
cause due to the angular momentm conservation the ro- 
tational perturbation always decays in expanding phase; 
this situation will be reversed in the collapsing phase 
where the rotational mode will grow again due to the 
angular momentm conservation. We will consider a flat 
background with two reasons. First, even in the linear 
perturbations it is known that the relativistic-Newtonian 
correspondence is rather ambiguous in the presence of the 
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background curvature; this is the case if we include the 
gravitational potential in making the correspondence [8] . 
Second, current observations of the large-scale structure 
and the cosmic microwave background radiation favour 
near flat Friedmann world model with non- vanishing cos- 
mological constant [9]. In fact, we will include the cos- 
mological constant in our analyses, thus our results are 
relevant to currently favoured cosmology. 

Our previous study on the second-order perturbations 
revealed which gauge condition (equivalently, gauge- 
invariant combinations) suits our problem [2]. In fact, 
unique gauge conditions were distinguished in showing 
the correspondence: these were the spatial C-gauge and 
the temporal comoving gauge. The spatial C-gauge was a 
rather natural choice because only in this gauge the spa- 
tial gauge mode can be fixed completely to all orders. In 
this gauge all the rest of the variables can be equivalently 
regarded as spatially gauge-invariant ones to all orders, 
see [1]. In the temporal gauge choice, however, we have 
many (in fact, infinitely many) different gauge choices, 
which fix the temporal gauge mode completely. In any 
of such gauge conditions all the rest of the variables can 
be cqiiivalently regarded as temporally gauge-invariant 
ones. Such procedures for taking the temporal gauge 
condition are separately available to each order in per- 
turbations; i.e., we can choose different gauge conditions 
to different perturbational order, see [1]. 

The (temporal) comoving gauge distinguished itself 
with the following reasons. In [2] we have successfully 
shown that in the comoving gauge we can identify the 
density and velocity variables which allow us to derive rcl- 
ativistic equations identical to the Newtonian ones. This 
point, perhaps, does not necessarily imply that it is not 
possible to discover any other variables and gauges which 
also lead to the same identification. But, we do mean that 
the comoving gauge is a natural choice. This also docs 
not necessarily imply that even to the third order the 
comoving gauge choice will be the best choice. As a mat- 
ter of fact, our policy /strategy about the gauge choice is 
that we do not know which gauge will suit the problem 
before we investigate and try many different gauge con- 
ditions. That is why we have presented our basic set of 
equations in a gauge-ready form which allows maximal 
usage of many different gauge conditions: see [10,11] for 
the linear case and [1] for the second-order perturbations. 

In this work, based on our successful experiences in 
the linear and second-order perturbations we will take 
the temporal comoving gauge condition and the spatial 
C-gauge even to the third order. First, under this gauge 
condition we can write the fully nonlinear equations in 
simple forms. Second, only in this gauge we have the 
relativistic equations the same as the Newtonian ones 
to the second order. And third, quite interestingly, we 
will show that in order to derive the third-order pertur- 
bations in this gauge, in fact, we need to evaluate the 
geometric and energy-momentum variables to the sec- 
ond order only. This is possible due to our appropriate 
choice of the variables and the gauge conditions. Thus, 



in the following we will derive the third-order pure rel- 
ativistic correction terms in the comoving gauge condi- 
tion. As we mentioned, since our gauge conditions fix 
the gauge modes completely, each of the variables we are 
using has a unique corresponding gauge- invariant combi- 
nation. Thus, our analyses can be equivalently regarded 
as gauge-invariant ones extended to the third-order per- 
turbations. 

Here, we summarize the relativistic-Ncwtonian corre- 
spondences known up to second-order perturbations in 
a flat Friedmann world model without pressure. To the 
background order we have [3,4] 
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with the energy (mass) density /i [g] oc a^^] a{t) is the 
cosmic scale factor and H = d/a; we set c = 1. The 
"const." part is interpreted as the spatial curvature in 
Einstein's gravity [3], and the total energy in the New- 
ton's gravity [4]. We put the cosmological constant A by 
hand which can work as a repulsive force proportional to 
the distance for A > 0. To the linear-order perturbations 
we have [5,6] 



5 + 2HS - inGgS = 0, 



(2) 



where S = S/i/fj, = 6q/q with fj, (g) and Sfi {5g) the back- 
ground and perturbed parts of the energy (mass) density. 
This equation is valid considering general presence of the 
background spatial curvature term and the cosmological 
constant. Now, to the second order, assuming flat back- 
ground, we have [7,1,2] 



(5-|--V-u = --V- (Su) 
a a 

1, 



(3) 



V ■{u + Hu)+4TrGgaS = --V-{u-Vu), (4) 
or by combining these, we have 



6 + 2H5 - A-kGqS 
1 d 



-^^[aV.(H] + ^V.(u.Vu). 



(5) 



These equations are valid in the presence of cosmological 
constant. The above equations are valid in both Ein- 
stein's and Newton's theories. In the relativistic theory 
we have to specify the variables 5 and u which correspond 
to the relative density fluctuation and perturbed velocity 
in Newton's theory; this will be done later. 

We believe no one would have anticipated such an ex- 
act coincidence to the second-order perturbations espe- 
cially considering the presence of the horizon in the rel- 
ativistic treatment. It might happen as well that our 
relativistic results give relativistic correction terms ap- 
pearing to the second order which could become impor- 
tant as the scale approaches and goes beyond the horizon. 
Our results show that there are no such correction terms 
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appearing to the second order, and the correspondence is 
exact to that order. Equations (3)- (5) are vahd in fully 
nonlinear situation in Newton's theory [7] whereas these 
are valid only up to the second order in Einstein's case 
[1,2]. It is our task to derive the third-order correction 
terms in Einstein's theory and to show how it causes dif- 
ference between the two theories. 



II. FULLY NONLINEAR EQUATIONS 

In [2] we have presented the fully nonlinear equations 
in the comoving gauge condition using the 1 + 3 covariant 
formulation [12] and the 3 + 1 ADM (Arnowitt-Deser- 
Misner) [13] formulations of Einstein's gravity. In [2] 
we showed that these two formulations are equivalent, 
and in the following we will take the ADM formulation. 
The basic set of ADM equations in our notation can be 
found in §11. A of [1], see also [14]. As this work can be 
regarded as a continuation of our previous studies in [1,2], 
we follow the notations used in that works. 

In the ADM approach, the temporal comoving gauge 
condition to all orders sets the flux four-vector to vanish, 
i.e., Ja = 0; here we also used the irrotational condition 
which ignores the vector-type perturbation. In [1] the 
fluid quantities are introduced based on the normal frame 
four-vector fia with fia = 0; in this case the information 
of the fluid motion is present in the flux vector ija with 
Qafi'^ = 0. In such a choice of the frame the temporal co- 
moving gauge condition with vanishing rotation implies 
Ua = of the fluid four- vector Ua- Thus, the fluid four- 
vector coincides with the normal four- vector. The phys- 
ical zero-pressure condition implies vanishing isotropic 
pressure (p or S in the ADM notation) and anisotropic 
stress {wab or ^a/j), i.e., p = = jtab or S = = Sa/s 
to all orders, based on the fluid four-vector. As the Mq 
coincides with the in our comoving gauge, the zero- 
pressure conditions simply allow us to set all the pressure 
terms in the normal-frame fluid quantites equal to zero; 
we point out that this is not true in the other gauge con- 
ditions in the normal- frame fluid quantities, see [2]; in the 
normal frame, in other than the comoving gauge condi- 
tions the pressure and anisotropic stress do not vanish 
to the second and higher orders even in the physically 
zero-pressure situation; for further discussions, see [2]. 
Although we do not need the form of energy-momentum 
tensor, from Eqs. (3), (4) of [1] we have 
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N-^E, f-o = 0, t 



afi 



0, 



(6) 



where tildes indicate the covariant quantities; E is the 
ADM energy density and N is the lapse function defined 
as 7V2 = -(gon)-i. 

The momentum conservation equation in Eq. (13) of 
[1] gives 



0. 



(7) 



Thus, we may set N = a(f) to all orders. In this case 
we have E = E^qN~^. The energy conservation equation 
and the trace part of ADM propagation equation in Eqs. 
(12), (10) of [1] give 



E-KE = Q, 



(8) 



k - -K'^ - K"''^K^i3 - A-kGE + a = 0, (9) 

o 

where E = E - E^aN^N^^, etc.: K and K^fi arc the 
trace and tracefree parts, respectively, of the extrinsic 
curvature of the normal hypersurface introduced in 
the ADM formulation, and Na is a shift vector defined as 
Na = (jQa- The spatial indices in the ADM formulation 
are based on the ADM three-space metric hajs defined as 
hoip = 5a/3- By combining these equations we have 



K^'^KafS - AttGE + a = 0. (10) 



Equations (6)-(10) are vahd to all orders, i.e., these equa- 
tions are fully nonlinear. 



III. THIRD ORDER PERTURBATIONS 

We consider the scalar- and tensor-type perturbations 
in the flat Priedmann background. As the metric we take 




-a^ (1 + la) dif - 2a^l3^ccd'ndx° 



+a^ 



5^ (1 + 2<p) + 27,„|^ + 2CWl dx'^dx^ (11) 



where a, /3, 7 and are spacetime dependent perturbed- 
order variables; we take Bardeen's metric convention in 

[10,11] extended to the third order. A vertical bar indi- 

cates a covariant derivative based on which can be 
regarded as 6(xis if we use Cartesian coordinates in the 
flat Friedmann background. By taking 7 = 0, which we 
call the spatial C-gauge, the spatial gauge mode is re- 
moved completely, thus all the remaining variables we 
are using are spatially gauge-invariant to the third order; 
this is true if we simultaneously take a temporal gauge 
which removes the temporal gauge mode completely, see 
§VI.B.2 and C.l of [1]. In the following we will take 7 = 
as the spatial gauge condition and use x = a/3+a^7 which 
becomes x ~ 
We expand 



E = fx + Sn, K= -3H + K. 



(12) 



Up to this point our notations look exactly the same as 
in the linear theory whereas, in fact, we consider each 
perturbation variable to be expanded to the third order. 
As an example, for S/j. we have 



(13) 
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where, to the third order we truncate the expansion at 
third-order term dfjMK A close examination of our fuUy 
nonhncar equations in §11 reveals an important technical 
methods of performing the third-order perturbations by 
using only the second-order expansion of perturbation 
variables. Such a simple method is possible due to our 
right choice of the equations, the right gauge conditions, 
and our proper choice of the fundamental variables to 
be matched with the Newtonian variables to the third 
order. That is, in our calculations we do not even need 
to have the inverse metric expanded to the third order, 
and all algebraic quantities we need can be found in [1] 
which presents various useful quantities expanded to the 
second order in perturbations. 

To the linear order we identified [8] 



u = — Vu 



1. 



X' 



— Vu = 



To the second order we identified [2] 



5iJby = dg, Ky = — V • u. 

a 



(14) 



(15) 



Based on our experience in the second-order perturba- 
tions, we have in mind to identify Eq. (15) even to the 
third order. It may turn out to be that these are not 
the best identifications, but in the following we will as- 
sume these arc the right ones and will take the consequent 
additional third-order terms as the pure relativistic cor- 
rections. 

The perturbed-order variable v is defined as Ja = 
—aixv^a to all orders, thus v is related to the velocity or 
flux variable. Our comoving gauge condition sets u = 0. 
In our notation, Sfiy indicates a gauge-invariant combi- 
nation which is the same as Sji in the comoving gauge 
which sets 7; = 0. Such a variable (equivalently, a gauge- 
invariant combination) is unique to all orders in pertur- 
bations; to the linear order we have Sfjby = S/j, — ftav, and 
to the second order see Eq. (282) in [1]. Similarly, (p^ is a 
gauge-invariant combination equivalent to ip in the zero- 
shear gauge which sets x = 0; to the linear order we have 
(fij^ = (fi — Hx, and to the second order, see Eq. (280) in 
[1]. For our justification to name the gauge conditions, 
see below Eq. (265) of [1]. 

One of the terms we need to evaluate to the third order 
in Eqs. (8)-(10) is K"^Ka0- Since Kajs is already at least 
linear order, we need to evaluate K"^ to the second order 
only, etc. The other ones we need to evaluate to the third 
order are 



E = E-E^N°'N- 



K = K - K r.N°'N- 



(16) 



Due to our identifications in Eq. (15), we are using the 
perturbed parts of E and K in Eq. (12) as the fundamen- 
tal perturbation variables. In order to evaluate E,aN"-, 
since E^a (or iV") is already at least linear order, it is 



enough to evaluate TV" (or E^a) only to the second or- 
der, and similarly for K^aN". Thus, using Eqs. (55), 
(57), and (175) of [1], we 'have 

E,aN^ = --<5/z,„x'" (1 - 2^) + 2-6fi-'-x'''A (17) 
a a 

K,o.N" = aX'" (1 - 2^) + 2-K'"x'''C^l (18) 



a/3 



X 



X I (y^XMP + C^S) (1 - 2a - 4^) - ^x.aV'./s 



Thus, Eqs. (8), (9) give 



(19) 



^ + 3i? (1 + 5) 



= k5- ^^,aX'" (1 - 2^) + 2^5'«x'''C'S, (20) 
-m - ZH'^ - A-kGijl + k + k + 2Hk - AirGSfi 
= - ^« (1 - 2^) + 2^^n"^x^C% 

^(x'"l^-^5(^^"^Ax)+C(*)"'^" 

^^X,«|/3 + C^S) (l-2a-4</.) 
~XM - ^^C% - Ix'-^ (2C% - C^iJ 



-4C(*)xf^X,m + cS 



(21) 



We note that a and in the comoving gauge are already 
quadratic order at least, thus can be ignored in Eq. (21); 
see Eqs. (12), (20) of [2]. 

Now, we need Newtonian expressions of ip to the linear 
order and x to the second order. Since we are considering 
the comoving gauge condition our <^ and x are the same 
as the gauge-invariant combinations ^py and Xvi respec- 
tively. To the linear order we have 



Lpy = ip — aHv = 99^ — aHv^, 
Xv=X- = -a''x- 



(22) 



To the second order, Xv is presented in Eq. (284) of [1]. 
To the second order, from Eq. (197) of [1] we have 



A 



3 1 



is) 



V a' 



(2<pAx - X'>,a) 
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a 



(23) 



where we have ignored a and ip which contribute to the 



third order. Apparently, we also need C^l to the second 



order. This will be presented in §V. 



IV. SCALAR-TYPE CORRECTIONS 

Ignoring C^*^, the perturbed parts of Eqs. (20), (21), 
and (23) give the complete set. As we ignore the rota- 
tional type perturbation the Newtonian velocity pertur- 
bation u is of a potential type, i.e., u = Vu. Thus, Eqs. 
(15), (23) give 

K= ■u= --Au, x = au + aA~'^X, (24) 
a a 

where k is valid to the third order and x is valid to the 
second order. Thus, Eqs. (20), (21), and (23) can be 
written as 

(5+ iv • u = --V • (Su) 
a a 



1 



-V ■ (u + Hu) + AnGfiS = -^V • (u • Vu) 



(25) 



■-^(^u-V(V-u) + 4v 
3a^ 



<^ ( u • Vu — ■ " 



-^[u-V(A-i^)]+^u.VX + Axv-u, (26) 
where 

3 

X = 2<^V • u - u ■ V<y9 + 2^"^^ ■ [u • V (Vip) + uA(p] . 

(27) 

Equations (25), (26) extend Eqs. (3), (4) to the third 
order. By combining Eqs. (25), (26) we can derive 

S + 2-d- 4nGn5 = — ^— [aV ■ (^u)] + • (u ■ Vu) 
a at 



a? dt 

2 .4 



+_<^u.V(V-u)--;^V 



( u • Vu — gUV • u 



+ ^['^-V(A-^^)]-^u.VX- Axv-u, (28) 

which extends Eq. (5) to the third order. The last three 
lines of Eq. (28) are the third-order terms. Examination 
of Eqs. (25)-(27) shows that all the third-order correction 



terms are of 95-order higher than the second-order terms 
like -V ■ {5u) and Jj-V • (u ■ Vu). As the fully nonlin- 
ear zero-pressure Newtonian equations are exact to the 
second order, the above third-order correction terms in 
our relativistic analyses are pure relativistic correction 
terms. Therefore, the pure general relativistic effects are 
at least lyj-order higher than the Newtonian ones. Up to 
the third-order corrections appearing in the general rel- 
ativity, the effects arc independent of the horizon scale 
and depend on the linear-order curvature (gravitational 
potential, see below) perturbation strength only. 

In Eqs. (25)-(27) we need only to the linear order. 
Thus, let us examine the behavior of to the linear order. 
In our comoving gauge condition (f is equivalent to a 
gauge-invariant combination ipv, and to the linear order 
from Eq. (22) we have 



(29) 



where we have ip^ = —5$ and u = — Vw^ in Eq. (14). 
Thus, in terms of the Newtonian variables we have 



(p = -6<^ + dA~^V ■ u. 



(30) 



Exact solutions of linear perturbation were presented in 
Tables of [15]. For i^T = we have [15] 



fv = C, 

/* dt 

= aA-^K^ = -C-^ (^1 + a'HH j' ^ , (31) 

where the lower bounds of integrations give decaying 
modes. For K = = A, using a = ait^/^, we have 
[15] 

'■Pv = C, 

= ait^/'A-^K, = -- f _ ht-'A . (32) 

Notice that to the linear order we have 



(fiv = C, 



(33) 



and (fiv has no decaying mode in expanding phase; this 
is true considering the presence of the cosmological con- 
stant. In fact, to the linear order satisfies [8] 



(34) 



See also Eq. (20) in the second reference of [2] . Ignoring 
the decaying mode, for A = 0, we have 



(35) 
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and the temperature anisotropy of cosmic microwave 
background radiation (CMB) gives [16,17] 



ST 



1 



(36) 



to the Unear order. This is a part of the Sachs- Wolfe 

effect in a flat background without the cosmological con- 
stant. The observations of CMB give 5T/T ~ 10"^ [18], 
thus 



5 X lo-^ 



(37) 



in the large-scale limit near horizon scale. Our tp is di- 
mensionless. 

We call (p the ciirvature perturbation because it is re- 
lated to the perturbed part of the spatial curvature of 
the normal hypersurface. To the linear order we have 



6K A + 3K 
—2 2 — 



(38) 



where R^^^ is the scalar-curvature of the three-space met- 
ric hap, and K is the sign of the background curvature, 
see Eq. (4) in [11]. To the second order, see Eq. (265) 
in [1]. In fact, from Eqs. (7), (55), and (175) of [1] 
we can easily check that (p characterizes the three-space 
Riemann curvature i?'^'*''"^^^ to all orders in perturba- 
tions; this is true assuming pure scalar-type perturba- 
tion. The gauge-invariant combination (fy is also known 
to be one of the large-scale (super-sound-horizon) limit 
conserved variables in various situations including time 
varying equation of state, field potential, and generalized 
gravity theories [19]. It is also known to be conserved in 
the large-scale (supcr-sound-horizon) limit even in non- 
linear situations [20,1]. 



V. INCLUDING THE GRAVITATIONAL WAVES 

In this section we present the complete set of equa- 
tions to the third order, now including the contribution 
of gravitational waves, and also the gravitational wave 
equation complete to the second order. Using Eqs. (14), 
(15), and (24), Eqs. (20), (21), and (23) give 

d+-V-u = --V- (Su) 

a a 



-- [2ipu - V (A-^X)] • V6 + -S'"u^A 



(39) 



-V • {u + Hu)+4TTGfi5 



3a 



2 4 
!/?u • V (V • u) 



( u • Vu — • u 



+4/^iy«l/3 + C7(*)«/5 



/37 



1 



3a ^"/'^ 
where 



(*) 

7a|/3 



(40) 



X = 2ipV ■u-n-WLp+ 2^"^^ ■ [u • V (V</9) + uAip] 



+ 



01- 



(41) 



These generalize Eqs. (25)- (27) in the presence of the 
gravitational waves. By combining Eqs. (39), (40) we 
have 

5 + 2H5 - AttG^jlS 



(*) 

'a0 



1 ^ 
2 

A 



-^,,u.V(V.u)-lv. 



( U • Vu — • U 



[u.V(A-ix)]-lu.VX- Axv-u 



(42) 



which generalizes Eq. (28) to include the gravitational 
waves. 

Now, we present the equation for tensor-type pertur- 
bation to the second order. From Eq. (103), (210) of [1] 

we can derive equation for C^*^ to the second order. Since 
we are ignoring the vector-type perturbation, from Eqs. 
(211), (199) of [1] we have 
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-I (v«V;3 - iff^A^ A-^W-^W'Ni^s. (43) 
From Eq. (103) of [1] to the second order we have 



\^ l(a.\p)o at}\~io 70|a/3 J 



a/317 



4%''''' X,7|(5 



kAx - 4^A<p - ^ip'^ip,^ + 2(/p'Tl'5c5 



_4C(*)7^ + C7(*)T^N (2C^*), - 3C(;) ^) ] | . (44) 

In Eq. (44) we have ignored a and terms which are 
already quadratic order in the comoving gauge. We have 
X = Xv, "i^ = "-pv, 1^ = i^v and = cfj^^ which are 
gauge-invariant combinations. Apparently, we need Xvi 
Kv and ipv to the hnear order. To that order, we have 
identified k„ = — • u and u = ^Vxi-- For Lp^ we have 
Eq. (29). Using these identifications we can express the 
scalar-type perturbation variables in Eq. (44) in terms of 
the Newtonian variables. 



VI. DISCUSSION 

We have derived the third-order perturbation equa- 
tions in the zero-pressure cosmological medium in Ein- 
stein's gravity. We have expressed the third-order terms 
using the Newtonian variables identified in the lower or- 
der perturbations. Since the Newtonian zero-pressure 
cosmological medium is exact to the second order in per- 
turbation, our third-order terms in relativistic context 
are pure relativistic corrections. Our results show that 
the third-order correction terms in relativistic energy and 



momentum conservation equations arc of ^s^-ordcr higher 
than the second-order terms, thus cquivalently (/^^-order 
higher than the Newtonian terms. The corrections terms 
are independent of the horizon scale and depend only 
on to the linear order which is the spatial curvature 
perturbation in the comoving gauge (hypcrsurfacc) , or 
~ i5$ which is the gravitational potential perturbation 
strength. The variable ipv is known to have conserved 
behavior and its amplitude in the large-scale (near hori- 
zon, say) is constrained by the low-level anisotropics of 
the CMB temperature, see Eq. (37). Therefore, our re- 
sult reinforces our previous conclusion in [2] that one can 
use the large-scale Newtonian numerical simulation more 
reliably even as the simulation scale approaches near (and 
goes beyond) the horizon. 

In this work we have assumed a single zero-pressure 
irrotational fluid in the flat cosmological background. 
Dropping any of these conditions could potentially lead 
to relativistic corrections. Due to these assumptions we 
cannot apply our results when the radiation components 
(including neutrino anisotropic stress) become important 
in high redshift epoch, and in the case when the baryon 
generated entropy leads to rotational perturbations in 
the small-scale clusters.* Extensions to include relativis- 
tic second-order perturbational effects of the pressure, 
the rotation, the non-flat background, and the multi- 
component situation will be investigated in future occa- 
sions. In this work we derived the relativistic correction 
terms appearing in the third order, and showed that these 
correction terms do not involve the horizon scale and are 
small in our observable patch. 

Besides theoretical and practical signiflcance we be- 
lieve our exact result to the second order and pure rela- 
tivistic corrections to the third order have historical value 
as well, because these have been unsolved issues since 
Lifshitz's original work in the linear regime in 1946 [5]. 
Before our present third order and the previous second 
order works there were different anticipations among re- 
searchers in the field that even in the second order the 
relativistic result might be different from the Newtonian 
ones: one common anticipation was that the general rel- 
ativistic effects might become important as the scale ap- 
proaches and goes beyond the horizon. Our results re- 
solved such an issue and showed that there exist no cor- 
rection terms to the second order in all scales. We also 
showed that pure relativistic correction terms appearing 
in third order do not depend on the horizon scale. It de- 
pends only on the strength of the dimensionless gauge- 
invariant curvature perturbation variable ip^ or the di- 
mensionless linear-order gravitational potential 6^. 

The post-Newtonian approximation [21] provides a 
complementary method to our perturbative approach in 
deriving the relativistic correction terms in the New- 



*We thank the anonymous referee for making this point. 
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tonian cosmology. The post-Newtonian approximation 
takes t'/c-expansion with v/c <^ 1, thus for near virial- 
ized systems we have GM/{Rc^) ~ v'^/c^ < 1. Thus, 
although such an approximation takes into account of 
nonlinear ity it is valid only far inside the horizon; as we 
approach the horizon GM/{Rc^) becomes unity. Con- 
sistency of the Newtonian cosmology with the Newto- 
nian limit of the post-Newtonian approximation was re- 
ported in [22]. Wc can show that the Newtonian cosmo- 
logical hydrodynamic equations naturally appear in the 
zeroth-order post-Newtonian approximation [23]. Re- 
cently, we derived the fully nonlinear first-order post- 
Newtonian correction terms, and showed that these cor- 
rection terms have typically GM/{Rc^) l<? ~ 10~^ 
order smaller than the Newtonian terms in the non- 
linearly clustered regions [23]. Being a complimentary 
approach to the post-Newtonian approximation (which 
provides fully non-linear equations) our equations valid 
to the third-order perturbations may have diverse ap- 
plications in the cosmological situations where the sys- 
tems have not reached fully nonlinear stage. Compared 
with the post-Newtonian approach, our perturbation ap- 
proach is applicable in all cosmological scale. 

Even in the small (far less than the horizon) scale the 
pure (third-order) relativistic correction terms could have 
important roles if the strength of linear-order is large 
enough. Our relativistic results are valid in the pertur- 
bative sense. Thus, if (p„ approaches near unity, higher- 
order perturbative terms could become important as well, 
and it is likely that our perturbative approach breaks 
down. Still, it would be interesting to investigate regimes 
where ip^ is moderately important so that we can study 
the roles of pure relativistic effects using our third-order 
correction terms. For example, whether such pure rela- 
tivistic correction terms could lead to an observationally 
distinguishable non-Gaussian signature [24] is an inter- 
esting issue which may deserve further attention. For 
such investigations, Eqs. (25)-(28) are the complete set 
for pure scalar- type perturbation, and Eqs. (39)- (44) pro- 
vide the complete set including the gravitational waves. 
As we consider a flat background the ordinary Fourier 
analysis can be used to study the mode-couplings as 
in the Newtonian case in [25]. Comparing the roles of 
pure relativistic third-order corrections with the inher- 
ent third-order perturbation effects in the Newtonian ap- 
proximation in Eqs. (25)- (28) will be an interesting step 
we can take. Applications to such cosmological situations 
are left for future studies. 
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